Testing standard and non-standard neutrino physics with cosmological data 
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Cosmological constraints on the sum of neutrino masses and on the effective number of neu- 
trino species in standard and non-standard scenarios are computed using the most recent available 
cosmological data. Our cosmological data sets include the measurement of the Baryonic Acoustic 
Oscillation (BAO) feature in the Data Release 9 CMASS sample of the Baryon Oscillation Spec- 
troscopic Survey (BOSS). We study in detail the different degeneracies among the parameters, as 
well as the impact of the different data sets used in the analyses. When considering bounds on the 
sum of the three active neutrino masses, the information in the BAO signal from galaxy clustering 
measurements is approximately equally powerful as the shape information from the matter power 
spectrum. The most stringent bound we find is ^ m„ < 0.32 eV at 95% CL. When non-standard 
neutrino scenarios with N e a massless or massive neutrino species are examined, power spectrum 
shape measurements provide slightly better bounds than the BAO signal only, due to the breaking 
of parameter degeneracies. 

Cosmic Microwave Background data from high multipoles from the South Pole Telescope turns out 
to be crucial for extracting the number of effective neutrino species. Recent BOSS data combined 
with CMB and Hubble Space Telescope measurements give N e s = 3.66jlo'liio 69 m the mass less 
neutrino scenario, and similar results are obtained in the massive case. The evidence for extra 
radiation N e a > 3 often claimed in the literature therefore remains at the 2a level when considering 
up-to-date cosmological data sets. Measurements from the Wilkinson Microwave Anisotropy Probe 
combined with a prior on the Hubble parameter from the Hubble Space Telescope are very powerful 
in constraining either the sum of the three active neutrino masses or the number of massless neutrino 
species. If the former two parameters are allowed to freely vary, however, the bounds from the 
combination of these two cosmological probes get worse by an order of magnitude. 

PACS numbers: 98.80.-k 95.85.Sz, 98.70.Vc, 98.80.Cq 
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I. INTRODUCTION 

Over decades, solar, atmospheric, reactor, and accel- 
erator neutrinos have provided compelling evidence for 
the existence of neutrino oscillations, implying non-zero 
neutrino masses. Oscillation experiments only provide 
bounds on the neutrino mass squared differences, and 
therefore information on the overall scale of the neutrino 
mass must come from other experiments. Cosmological 
data provides a tool to test the absolute scale of neu- 
trino masses. Neutrino masses and abundances leave 
crucial features in several cosmological observables. The 
amount of primordial relativistic neutrinos affects Cos- 
mic Microwave Background (CMB) anisotropies and non- 
relativistic neutrinos in the recent Universe suppress the 
growth of matter density fluctuations and galaxy clus- 
tering. Cosmology can therefore weigh neutrinos, pro- 
viding an upper bound on the sum of the three active 
neutrino masses, m v ~ 0.26 eV at 95% confidence 
level (CL) [l| (see also Refs.@-Q). The former bound is 
found when CMB measurements from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) are combined with 
information on the distribution of galaxies based on a 
catalog of photometric galaxy redshifts in Sloan Digital 
Sky Survey 111 Data Release Eight (SDSS DR8), and with 
the most recent measurement of the Hubble parameter 
from the Hubble Space Telescope (HST). ft assumes a 
flat universe with a cosmological constant, i.e. a ACDM 



cosmology. 

In the Standard Model of elementary particles, there 
are three active neutrinos. However, there is no funda- 
mental symmetry in nature establishing a definite num- 
ber of right-handed (sterile) neutrino species. In fact, 
models with additional ~ 1 eV massive sterile neutri- 
nos [5j have been introduced to explain short baseline 
(SBL) antineutrino data [|| in terms of neutrino oscil- 
lations. Up to date cosmological constraints on massive 
sterile and active neutrino species have been presented in 
different cosmological scenarios using different data sets, 
see Refs. [Ml- 

In addition, extra relativistic degrees of freedom could 
also arise from different physics, such as massless ster- 
ile neutrino species, axions 1131. decay of non-relativistic 
matter [TJ], gravity waves [la], extra dimensions [llij . 
early dark energy [17| or asymmetric dark matter mod- 
els [10| . Cosmological probes have been extensively ex- 
ploited in the literature to set bounds on the relativis- 
tic energy density of the universe in terms of the ef- 
fective number of neutrinos N c g (see Refs. 0-0, Il9l - 
I22I]). Measurements of Cosmic Microwave Background 
(CMB) anisotropies at arc-minute angular scales from 
the South Pole Telescope (SPT) [H] and the Atacama 
Cosmology Telescope (ACT) [2J], when combined with 
other cosmological data sets have provided the constraint 
N cff = 4.08±^g at 95% CL [HHl), showing evidence 
for N e g > at more than la and suggesting values higher 
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than those expected in the canonical scenario. If the ef- 
fective number of neutrino species N c g is larger than the 
Standard Model prediction, i.e. if N c g > 3.046 during the 
Big Bang Nucleosynthesis (BBN) era, the Hubble expan- 
sion rate will be larger, causing therefore weak interac- 
tions to freeze out earlier. The standard BBN predictions 
for light element abundances will change. Current analy- 
ses of the observed primordial abundances indicate best- 
fit values of AN cS ~ 0.5 - 0.8 @, ©, H, HI, and bounds 
AN e g < 1 — 2 at 95% CL, depending on the assumptions 
on the fiducial cosmology, and are consistent with the 
canonical AA^ff — 0. Nevertheless, the extra radiation 
species do not necessarily feature thermal abundances 
at BBN [3l], HH so that the true number of additional 
species could possibly be larger than the effective num- 
ber defined in terms of thermal abundances. 

Here we follow a general approach and study the con- 
straints on several neutrino scenarios with the most re- 
cent available cosmological data. We consider both stan- 
dard neutrino scenarios with three active massive neutri- 
nos as well as more general, non-standard schemes with 
N e g massless or N c g massive species. Of the two sce- 
narios with free number of neutrino species, N e g, the 
latter one is the most realistic and appropriate scenario, 
since we know observationally that neutrinos have mass. 
However, this scenario has been explored much less in the 
literature than the case of N e g neutrinos with mass fixed 
to zero (although see, e.g., @, and the very recent 
[l2l]). In fact, we will show that in the scenario with N c g 
massive neutrino species, the neutrino constraints from 
some combinations of cosmological data sets are much 
worse than in the other two simpler schemes where one 
of the two parameters is kept fixed. This highlights the 
danger of neglecting neutrino mass. 

In each of the three scenarios analysed here, we shall 
study the impact of each data set separately, devoting 
special attention to parameter degeneracies and to the 
impact of galaxy clustering data. We shall also study the 
preferred values of N Q g arising from our different numeri- 
cal studies. As stated above, N e g has been claimed to be 
larger than the standard model expectation in a plethora 
of cosmological data analyses. We shall see that using up- 
to-date cosmological data, the evidence for N e g > 3 still 
remains at the 2a level. The data sets exploited include 
recent measurements of the BAO feature in the Data Re- 
lease 9 (DR9; IH) CMASS sample of the Baryon Oscilla- 
tion Spectro scop ic Survey (BOSS) [34],[35|], which is par t 
of SDSS-III [if with a median redshift of z = 0.57 

The structure of the paper is as follows. Section |TT] 
describes the three standard and non-standard physics 
scenarios considered here, summarizing the main effects 
of the neutrino parameters on the cosmological observ- 
ables. In Sec. IIIII we describe the data sets used in the 
numerical analyses and in Sec. HVI the constraints on the 
neutrino thermal abundance N c g and on the sum of their 
masses are presented for different combinations of data in 
the three schemes studied here. We draw our conclusions 
in Sec. EJ 



II. THE MODEL: NEUTRINOS IN 
COSMOLOGY 

We shall start with the standard scenario, which is 
made of three massive neutrinos and put constraints 
on the sum of their masses ^2 m„ . We therefore fix 
N e g — 3.04. Current neutrino oscillation data can not 
distinguish between the two possible types of ordering of 
the neutrino mass eigenstates, i.e. information regarding 
the neutrino mass hierarchy is inaccessible. While the ef- 
fect of neutrino mass on the CMB is related to the total 
amount of energy density in the form of massive neutrinos 
and the individual values of the neutrino masses are irrel- 
evant, the matter power spectrum in principle provides 
information on the sum of the neutrino masses as well as 
on the individual values of the m Vi 's. The non-relativistic 
neutrino overdensities cluster at a given redshift z only at 
scales where the wavenumber of perturbations is below 
the neutrino free streaming scale 

<=/.w = (i^(fw)^" M p«-'. m 

due to the pressure gradient, which prevents gravita- 
tional clustering (Q m is t ne ratio of the total matter en- 
ergy density over the critical density at redshift zero). 
On spatial scales larger than the free streaming scale 
k < kf g , neutrinos cluster and behave as pressureless 
matter (dust). Perturbations with comoving wavenum- 
ber larger than the free streaming scale can not grow 
due to the large neutrino velocity dispersion. As a con- 
sequence, the growth rate of density perturbations de- 
creases and the matter power spectrum is suppressed at 
k > kt 3 . Therefore, there is a clear signal on the mat- 
ter power spectrum induced by the presence of neutrino 
masses. Neutrinos of different masses will have different 
transition redshifts from relativistic to non-relativistic 
behavior, and, in principle, by exploiting the information 
contained on the matter power spectrum it could be pos- 
sible to identify the neutrino mass hierarchy and isolate 
the individual neutrino masses. However, in practice, 
such a task has been shown to be extremely challeng- 
ing [381 ] . In the following, we will assume a degenerate 
mass spectrum. In the standard scenario with three ac- 
tive massive neutrinos, the parameters considered in the 
analysis are: 

{w 6 ,u c ,e„T,n„log[10 10 i4a],5^m v } , (2) 

where w& = flbh 2 and uj c = fl c h 2 are the physical baryon 
and cold dark matter densities, O s is the ratio between 
the sound horizon and the angular diameter distance at 
decoupling, r is the optical depth, n s is the scalar spectral 
index, A s is the amplitude of the primordial spectrum 
and ^2 vn v is the sum of the active neutrino masses. 

Next, we explore non-standard neutrino scenarios in 
which the effective number of thermalized species is pa- 
rameterized by A/ejf. As commonly assumed in the lit- 
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erature, we start here assuming a massless neutrino sce- 
nario with N e ff massless species. One of the main effects 
of N e s on the CMB temperature anisotropies arises from 
the change of the epoch of the radiation matter equal- 
ity, shifting therefore the location of the CMB acoustic 
peaks. This position is given by the so-called acoustic 
scale 8a, which reads 

e A = r -^r4 , (3) 

where rg(z rec ) and r s (z rec ) are the comoving angular di- 
ameter distance to the last scattering surface and the 
sound horizon at the recombination epoch z rec , respec- 
tively. Although rg(z rec ) almost remains the same for dif- 
ferent values of N c fi, r s (z rec ) becomes smaller when JV e ff 
is increased. Therefore, the positions of acoustic peaks 
are shifted to higher multipoles (smaller angular scales) 
if the value of N c g is increased. However, this effect can 
be compensated by changing the cold dark matter den- 
sity, in such a way that z rec remains fixed, see Ref. [26| . 
Therefore, due to the degeneracy with the cold dark mat- 
ter component, the change induced at low I is negligible 
and the largest impact of N e g on the CMB temperature 
anisotropies comes from its effect on high multipoles £, 
since a higher value of N c g will induce a drop in power 
at small scales due to an increased Silk damping. Silk 
damping refers to the suppression in power of the CMB 
temperature anisotropies on scales smaller than the pho- 
ton diffusion length. As a random walk process, the dif- 
fusion distance r^ will increase as the square root of time 
and therefore a higher expansion rate H(z) (caused by 
a higher -/V e ff) will decrease r<j, decreasing therefore the 
damping. However, in order to keep the acoustic scale 
9a fixed, the comoving angular rg(z rec ) distance at de- 
coupling should scale in the same way as r s (z rec ) does 
(i.e. as 1/H(z)) and the angular scale of the Silk damp- 
ing 9d = rd/rg(z rec ) will grow as y/~H. Consequently, 
the damping will be increased when N^n does due to the 
higher expansion rate H(z), see Ref. 26]. 

This first non-standard scenario is described by the 
following set of parameters: 



with equal 1 mass m v . The effects of these two parameters 
on the observables have been summarized above. While 
the highest sensitivity to the neutrino masses comes from 
large scale structure information, and the highest N e g 
sensitivity comes from CMB measurements, there will be 
significant covariance between the two parameters and 
we explore this degeneracy in IIVB 21 In this last non- 
standard scenario, the set of parameters involved in the 
study are 

For our numerical analyses, we have used the Boltzmann 
CAMB code [4(| and extracted cosmological parameters 
from current data using a Monte Carlo Markov Chain 
(MCMC) analysis based on the publicly available MCMC 
package cosmomc [4l[. Table Q] specifies the priors consid- 
ered on the different cosmological parameters 2 . Our neu- 
trino mass prior is cast in the form of a (uniform) prior 
on the neutrino density fraction /„ = fi y /^DM ; wnere 
il„ is the ratio of the neutrino energy density over the 
critical density at redshift zero, and £2dm is the same ra- 
tio, but for the total dark matter density, which includes 
cold dark matter and neutrinos. 

We restrict ourselves to a ACDM cosmology in the 
different scenarios described above. 



Parameter 


Prior 




0.005 ->■ 0.1 


n c h 2 


0.01 -»■ 0.99 




0.5 -> 10 


T 


0.01 0.8 


n s 


0.5 -)■ 1.5 


ln(10 10 7l s ) 


2.7 ->■ 4 


u 


-> 0.2 


AW 


1.047 10 



TABLE I: Uniform priors for the cosmological parameters 
considered here. 



III. DATA 



{uj b , uj c , 6 S , r, n„ logllO 10 ^], N eS } 



(4) 



However, we know from neutrino oscillation experi- 
ments that neutrino have non-zero masses and mixings. 
In fact, a recent neutrino oscillation analysis has shown 
that models with three active plus two sterile light (sub- 
eV) neutrino species provide a very good fit to short 
baseline data [30|. Therefore the previous minimal non- 
standard scenario with N e g massless neutrinos should 
be enlarged to accommodate active plus sterile neutrino 
masses and mixings. The last scenario explored here 
therefore consists of N e g massive neutrino species, each 



Our baseline data set is the seven-year WMAP data 
2, 42] (temperature and polarization) with the routine 
for computing the likelihood supplied by the WMAP 
team. We then also add CMB data from the South Pole 



1 Also in the presence of additional species, current data are not 
sensitive to how the masses are distributed over the species, but 
only to the sum of the masses and the number of species. This 
justifies the use of a degenerate mass spectrum. 

2 In the next section we will describe the additional nuisance pa- 
rameters that we have also added in the analyses. These param- 
eters are related to some cosmological data sets exploited here. 
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Telescope (SPT) [23|, which strongly improve the mea- 
surement of the temperature anisotropies on scales < 10 
arcmin. In order to address foreground contributions, the 
SZ amplitude A$z, the amplitude of the clustered point 
source contribution, Ac, and the amplitude of the Pois- 
son distributed point source contribution, Ap, are added 
as nuisance parameters in the CMB data analyses. We 
have followed Ref. [23j , applying a Gaussian prior on the 
amplitude of each of these three foreground terms. 

To the WMAP basic data set we add the latest con- 
straint on the Hubble constant Hq from the Hubble Space 
Telescope (HST) [43[ , and supernova data from the 3 year 
Supernova Legacy Survey (SNLS3), see Ref. [44|. In the 



case of SNLS3 data, we add in the MCMC analysis two 
extra nuisance parameters related to the light curve fit- 
ting procedure used to analyze the supernova (SN) data. 
These parameters characterize the dependence of the in- 
trinsic supernova magnitude on stretch (which measures 
the shape of the SN light curve) and color [44[ . 

We also consider data from galaxy clustering that we 
summarize next. We first employ data from SDSS-II 
[Sloan Digital Sky Survey; [Hf) Data Release 7 (DR7; 
146]), analyzing separately the full shape of the halo 
power spectrum derived from the clustering of luminous 
red galaxies [47|, and the Baryon Acoustic Oscillation 
(BAO) signal as extracted from the same data set [48j . 
Then we consider the impact of the most recent mea- 
surement of the BAO scal e 1371 from the CMASS sample 
in Data Release 9 (DR9; [3j]) of the Baryon Oscillation 
Spectroscopic Survey (BOSS) |34|,l35j], with a median red- 
shift of z = 0.57. Together with the CMASS DR9 data, 
we also include the recent measurement of the BAO scale 
based on a re-analysis (using reconstruction |49|) of the 
LRG sample from Data Release 7 with a median redshift 
of z = 0.35 [Hoj], and the measurement of the BAO sig- 
nal at a lower redshift z = 0.106 from the 6dF Galaxy 
Survey 6dFGS [Hj]. We will refer to the combination of 
these three new BAO signals as BAO2012- 



IV. RESULTS 

Here we present the constraints from current data on 
the neutrino thermal abundance N c g and on the sum of 
their masses in different scenarios. 



A. The standard lore: Three massive neutrinos 

In this section we consider the standard picture, 
i.e. three massive neutrinos and put constraints on the 
sum of their masses m v . 

Table [H shows the 68% and 95% CL bounds on £ m v 
from several combinations of WMAP data with SPT, 
HST and SNLS3. These constraints agree quite well with 
previous analyses in the literature [l[ . Note that WMAP 
alone sets a 95% CL bound of ~ 1 eV on the sum of neu- 
trino masses. When HST data is added in the analysis 



the results improve in a very significant way, since HST 
data helps enormously in breaking the strong degener- 
acy between ^m u and Hq, see Fig. [TJ Indeed, with- 
out a Hq measurement, the change induced in the CMB 
temperature anisotropics caused by an increase in ^2 m v 
can be compensated by a decrease in Hq. The reason 
is the following: if ^2 Tn v increases, the main effect on 
the CMB comes from the resulting shift in the distance 
to last scattering 3 . While the acoustic peak structure of 
the CMB data does not leave much freedom in uj c and u>b, 
the change in distance can be compensated by lowering 
h (i.e. Hq). The HST prior on the Hubble parameter will 
break this strong degeneracy, setting a ~ 0.4 eV 95% CL 
bound on the sum of the three active neutrino masses. A 
similar result but to a lesser extent arises when combining 
WMAP and SNLS3 data sets. Note that we follow here a 
conservative approach avoiding the combination of HST 
and SNLS3 data, since these two data sets are not totally 
independent. The addition of SPT data does not help in 
constraining "m v '■ taking into account the marginaliza- 
tion over the additional nuisance parameters, the high £ 
multipole data are not sensitive to neutrino mass, nor 
do they help break degeneracies to improve the neutrino 
bound and therefore we will not consider SPT data set 
in the rest of the analyses presented in this section. Of 
course this will not be the case when the effective num- 
ber of neutrino species is added as a free parameter in the 
analyses, since, in that case, SPT high £ measurements 
will become crucial, as the combination of WMAP and 
HST will not be enough to break the strong degeneracy 
between N e g and m v ■ 

Table Hill shows the constraints in the standard massive 
neutrino scenario when SDSS-II galaxy clustering data 
are added in the analysis. Note that the bounds from 
geometrical information from the BAO signal are very 
similar to those obtained from the shape of the galaxy 
power spectrum measurements, although BAO informa- 
tion provides slightly better constraints on the sum of 
the three active neutrino masses. These findings agree 
with previous analyses in the literature, see e.g. Ref. 
However, we are considering here a minimal cosmological 
model, with spatial flatness and a cosmological constant. 
B AO / geometrical data will not be as powerful as shape 
measurements when other cosmological parameters are 
added in the analysis, since the information contained in 
the full shape of the matter power spectrum is extremely 
useful for breaking degeneracies. We will illustrate this 
effect in section HV B 21 where the addition of the num- 
ber of massive neutrino species N e ff will open up new 
degeneracies, and the matter power spectrum shape in- 
formation will be more powerful in resolving them than 
the geometrical BAO measurements. Note from Tab. IIIII 
that if other data sets are included, such as, for instance, 



3 rg(z ICC ) oc f*' 



dz [u> r a 4 +LO m a 3 + (1 — uim/h 2 
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L m 

V 

FIG. 1: The red contours show the 68% and 95% CL con- 
straints in the m„, Ho) plane from our basic WMAP data 
set. The blue contours show the results from the combination 
of WMAP and HST data, while the green contours depict 
the results from the combination of WMAP and SNLS3 data 
sets. Notice that the strong degeneracy present in the case of 
WMAP data alone gets alleviated when a prior on Ho from 
HST data is added in the analysis. SNIa luminosity distante 
data are unable to independently determine the Hubble con- 
stant Ho but measure the fl m quantity, and, since WMAP 
measures Q m h 2 , the combination of WMAP plus SNIa data 
is able to determine Ho- 




0.2^ J 

0.2 0.4 0.6 
L m 

V 

FIG. 2: The red contours show the 68% and 95% CL con- 
straints in the (J^m„, Q m ) plane from WMAP plus HST 
data sets. The blue contours show the results from the com- 
bination of WMAP, HST and galaxy clustering measurements 
from SDSS-II interpreted in the form of BAO signals, while 
the green contours depict the results from the combination of 
WMAP, HST and BAO2012 data sets. The errors on the O m 
parameter are significantly improved when considering either 
the SDSS-II BAO or the BAO2012 measurements. Also, due 
to the higher mean value of Q, m when adding BAO informa- 
tion to the WMAP and HST measurements, a slightly higher 
neutrino mass is allowed. 



a prior on Hq from the HST experiment, some degenera- 
cies are broken, and the constraints on ^2 fn u from shape 
and geometrical SDSS-II measurements of the galaxy dis- 
tribution become very similar again. In fact, once the 
HST prior is included, further inclusion of the large scale 
structure data considered here does not significantly im- 
prove the neutrino mass bound in the first place. 

Table IIVI shows the constraints on m„ with the ad- 
dition of the new BAO data from the CMASS DR9 BOSS 
experiment combined with a re-analysis of the LRGDR7 
sample and the BAO signal from 6dFGS. Note that the 
bounds in this minimal cosmological model from the 
BAO2012 data set combination are slightly better, albeit 
very close to those obtained with BAO SDSS-II measure- 
ments. Note from Tabs. Mil and IIVI that the combina- 
tion of WMAP and HST data with old BAO measure- 
ments gives a bound on ^2 m„ which is worse than the 
one arising from the combination of WMAP plus HST 
alone. The errors on the Cl m parameter arc significantly 
improved when considering either the SDSS-II BAO or 
the BAO2012 measurements, see Fig. [2] The mean value 
of f2 m is also higher when adding BAO information to 
the WMAP and HST measurements and, consequently, 
a slightly higher neutrino mass is allowed. 



B. Beyond the standard lore: varying N e g 

1, Neglecting neutrino mass 

We assume here a scenario with N e g massless neutrino 
species. Table IVl shows the central values and the 68% 
and 95% CL errors on the quantity N g for our basic 
WMAP data set alone and also with other data sets ex- 
plored here. We observe that WMAP data alone is un- 
able to constrain the value of N e g but the situation dras- 
tically improves when adding either SPT or HST data. 
The reason for that is because N e g is highly degenerate 
with both Hq and the dark matter energy density fluMh 2 
at the multipole range relevant to WMAP CMB mea- 
surements: the change induced in the CMB spectrum by 
iVeff can be compensated by a change in the value of the 
Hubble constant Hq. HST data is crucial for breaking 
the strong degeneracies between N s and both the phys- 
ical amount of matter f2oM^ 2 and the Hubble constant 
Hq, therefore, the prior on Hq from HST improves signifi- 
cantly the constraints on N e g. In the following section we 
will see when adding ^2 mv as a free parameter the errors 
on N c s from WMAP plus HST will be much higher, due 
to the strong degeneracy between the effective number of 
neutrino species and the sum of the neutrino masses. 

SPT data is also very powerful in constraining N e g 
due to the information contained at high multipoles I. 
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WMAP WMAP+HST WMAP+SPT WMAP+SNLS3 



£ m„ (eV) 68% CL <o.60 < 0.16 < 0.52 < 0.24 

95% CL < 1.05 < 0.37 < 1.07 < 0.51 



TABLE II: Constraints on m v from WMAP data alone and WMAP combined with HST measurements of the Hubble constant, 
SPT data and SNLS3 data. We have fixed the number of massive neutrinos to 3.046. 





WMAP+MPK WMAP+BAO WMAP+MPK+HST WMAP+BAO- 


1-HST WMAP+MPK+SNLS3 WMAP 


+BAO+SNLS3 


E m„ (eV) 68% CL 
95% CL 


< 0.37 < 0.29 < 0.16 < 0.22 

< 0.63 < 0.52 < 0.32 < 0.42 


< 0.22 

< 0.42 


< 0.27 

< 0.47 



TABLE III: Constraints on the sum neutrino masses m v from WMAP data and different possible combinations of galaxy 
clustering measurements from SDSS-II (interpreted either in the form of matter power spectrum measurements (MPK) or in 
the form of BAO signals (BAO)), HST and SNLS3 data sets. We have fixed the number of massive neutrinos to 3.046. 



As summarized in Sec. UH an increasing value of N e g 
has little (big) impact on low (high) multipoles I due to 
an increased Silk damping caused by a higher expansion 
rate. Therefore, SPT high multipole data helps enor- 
mously in constraining the value of N e s- The combina- 
tion of WMAP, SPT and HST gives 3.83±g |^±g;f | which 
deviates ~ 2a from the expected standard model value 
for N e ff. The combination of CMB measurements with 
SNLS3 luminosity distance data results in a mean value 
of N e g which is only ~ la away from its the standard 
value. 

The addition of the information contained in the mat- 
ter power spectrum or in the BAO signal from SDSS-II 
galaxy clustering measurements to our basic WMAP data 
set does not change significantly the central values nor 
the errors depicted in Tab. [V] and for the sake of clarity 
we do not show all the possible combinations here. We 
only illustrate the case with the addition of the new BAO 
data from the BOSS DR9 CMASS sample combined with 
a re-analysis of the LRGDR7 sample and the BAO signal 
from 6dFGS, see Tab. I VII The evidence for extra radi- 
ation N e g > 3 claimed in the literature still remains at 
the 2a level, in perfect agreement with a recent analysis 
presented in Ref. [l2| . 



2. Including neutrino mass 

We consider here a scenario with N e g massive species, 
each with a mass m v . Table IVIII shows the joint con- 
straints on the sum of neutrino masses ^2 m v an d on the 
number of massive neutrino species N c g when HST, SPT 
and SNIa data are added to our basic WMAP data set. 
Note that the addition of HST data to the WMAP basic 
data set does not help as in the massless neutrino sce- 
nario, the reason for that being that ra v is now a free 
parameter and it is strongly degenerate with Hq (among 
other parameters, such as the matter energy density). 
Therefore, a prior on Hq is not enough to improve the 



constraints from WMAP alone. The errors on N e g are 
an order of magnitude larger than in the massless neu- 
trino scenario. A similar effect occurs when addressing 
the neutrino mass bound: in the case where only three 
active massive neutrinos were considered, the WMAP 
and HST data combination provided a 95% CL bound 
m v < 0.38 eV. The former bound becomes 2.25 eV 
when iVcff is allowed to vary as well. 

High multipole SPT data helps in constraining N e g. 
SNIa data is very useful with respect to neutrino masses. 
However, both the mean value of N c g and its er- 
ror are very large (due to the strong degeneracy be- 
tween N e fi and ^2 m v as we shall see below) unless 
SPT data is added in the analysis. Such an addition 
(i.e. WMAP+SPT+SNLS3) allows to partially break the 
strong degeneracy between ^ m„ and N e g and it de- 
creases both the mean value of N e g and its errors sig- 
nificantly. Note from Table IVIII that the neutrino mass 
bounds from CMB measurements plus SNIa data are bet- 
ter than those from CMB plus HST measurements. After 
the addition of SNLS3 measurements, the data combina- 
tion prefers a lower value of J7 m which, in turn, leads to 
a smaller value of ^2 m v ■ 

Table IVIIII shows the constraints in the massive neu- 
trino scenario when data from galaxy clustering mea- 
surements are added in the analysis. When considering 
SDSS-II measurements, the bounds are stronger (about 
a factor of two) if the information from galaxy cluster- 
ing is exploited as shape constraints (i.e. power spectrum 
measurements) rather than as geometrical measurements 
(i.e. BAO signals), the reason for that being the presence 
of more degeneracies involved, as previously explained. 
When other data sets are included, some degeneracies are 
broken, and the constraints from shape and geometrical 
SDSS-II measurements of the galaxy distribution are very 
similar. We also depict the bounds on the sum of neu- 
trino masses and on the number of massive species when 
new BAO information from the BOSS DR9 CMASS sam- 
ple, combined with a re-analysis of the LRGDR7 sample 
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W M Ar+ D AU2012 


W M J\r-\-D AU2012 - 


(-rial WMAF+r$AU20l2+t3lNLibo 


£ m v (eV) 68% CL 


< 0.26 


< 0.19 


< 0.25 


95% CL 


< 0.48 


< 0.38 


< 0.47 



TABLE IV: Constraints on £ m„ from WMAP and different possible combinations of the new BAO data explored here 
(BAO2012), HST and SNLS3 data sets. We have fixed the number of massive neutrinos to 3.046. 



WMAP WMAP+HST WMAP+SPT WMAP+SPT+HST WMAP+SPT+SNLS3 



N eB 5.99 



+4.01+4.01 
-1.23-3.47 



4.19 



+0.33+1.30 



4.13 



-0.35H 
-0.41- 



1.46 
1.30 



3.83 



+0.21+0.86 



3.51 



+0.31+1.25 



TABLE V: Constraints on iV e fr for the massless neutrino scenario from WMAP data, WMAP data and HST measurements of 
the Hubble constant, WMAP and SPT data and WMAP plus SPT plus SNLS3 data. 



and with the BAO signal from 6dFGS, are added in the 
analysis. Note that these new BAO data sets are more 
powerful than the previous BAO data set because they 
combines BAO measurements at several redshifts, rely 
on a larger volume, and include an improvement of the 
original SDSS-II BAO measurement. Because of this, the 
BAO2012 constraints are even better than those obtained 
with the full matter power spectrum of SDSS-II. 

The strong degeneracy between N e g and £ m v is il- 
lustrated in Fig. [3l where we show the combination of 
WMAP with new and old BAO data sets as well as the 
combination of CMB data with HST or SNLS3 measure- 
ments. Finally, we also show the combination of CMB, 
BAO (both new and old data sets) and SNLS measure- 
ments. As in the massless neutrino case, the addition 
of new BAO2012 data to CMB and HST measurements 
makes the value of N e g (marginally) consistent with the 
standard expectation within 2a, in agreement with pre- 
vious analyses in the literature. Therefore the evidence 
for N e g > 3 still persists when considering massive neu- 
trino species. The 95% CL constraints on the effec- 
tive number of massive neutrino species and on the sum 
of their masses from CMB measurements, SNLS3 and 
BAO2012 data are 3.44±^ and ^ra„ < 0.47 eV, re- 
spectively. The former bounds translate into 3.71^o 72 
and J2 m v < 0.51 eV when CMB, HST and BAO2012 
measurements are considered. The bounds and the errors 
found in this work on both £ m v and N e g are similar to 
those presented in a recent study presented in Ref. [T^ |. 
where a preference for N c g > 3 at a 2a level has been 
presented. 



V. CONCLUSIONS 

We have presented constraints on several neutrino- 
extra radiation scenarios with the most recent and avail- 
able cosmological data. The data sets explored here in- 



clude the recent measurement of the BAO feature in the 
DR9 CMASS sample of the BOSS experiment. Three 
scenarios have been carefully examined: one standard 
scenario with three massive neutrino species, and two 
non-standard scenarios with N e g massless/massive neu- 
trino species. In each scenario we have carefully explored 
the different existing degeneracies and the impact of the 
different data sets exploited here, both separately and 
combined. In the standard scenario where only three ac- 
tive neutrinos are considered, the combination of CMB 
WMAP plus HST data helps enormously in breaking the 
strong degeneracy between £ m„ and the Hubble pa- 
rameter Hq that exists with CMB data only, and the 
addition of further cosmological data sets does not im- 
prove significantly the limit on £ m„ . In this standard 
scenario, WMAP data plus geometrical information from 
BAO signals gives similar results to those obtained from 
WMAP data plus the measurements of the full shape of 
the matter power spectrum. 

In the non-standard scenario with N e g massless neu- 
trinos, the combination of WMAP and HST data is also 
able to break the strong degeneracies between the effec- 
tive number of neutrino species and other cosmological 
parameters, such as the Hubble parameter and the physi- 
cal dark matter energy density SPT high multipolc CMB 
data is also extremely useful in constraining the effective 
number of neutrino species, due to the increase in the 
Silk damping effect at small scales induced in the case 
of iV c ff > 3. The addition of SDSS-II galaxy cluster- 
ing data or recent BAO measurements does not improve 
the constraints on N c g in the massless neutrino scenario. 
We obtain N cS = 3.66±g |5±g ^| when combining the for- 
mer BAO data set with CMB and HST measurements. 
Therefore, the cosmological evidence for N c g > 3 claimed 
often in the literature gets still remains when considering 
these new data sets, and the value of N c g very close to 
the standard model expectations (N e g = 3.04), being the 
discrepancy at the 2a level. 

In the most general non-standard scenario explored 
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WMAP+BAO2012 WMAP+SPT WMAP+SPT WMAP+SPT 

+BAO2012 +BAO2012+HST +BAO2012+SNLS3 



\T K QC+l-07+3.51 q ^cj+0.31+1.27 n fiR +0.20+0.73 o ro+< 

iV c ff °-°°-l. 15-3.36 J -^°-0.36-1.17 J - uu -0.21-0.69 



TABLE VI: Constraints on N c g from WMAP and different possible combinations of the new BAO data explored here (BAO2012), 
HST and SNLS3 data sets. 



WMAP WMAP+HST WMAP+SPT WMAP+SPT+HST WMAP+SNLS3 WMAP+SPT+SNLS3 



A7 R 1 0+3. 82+3. 82 r Q1 +4.09+4.09 o Q1 +0.30+1.32 , 1Q +0.25+1.02 c Q7 +4.03+4.03 n ,, R +0.32+ 

iv cft D - 10 -1.18-3.34 J - al -0.94-2.44 J ' 3± -0.39-1.15 ±3 -0.30-0.92 °' a ' -1.17-3.28 J '^ u -0.37- 

£m„ (eV) 68% CL < 0.90 < 1.07 < 0.81 < 0.46 < 0.28 < 0.22 

95% CL < 2.02 < 2.25 < 1.47 < 0.90 < 0.66 < 0.50 



TABLE VII: Constraints on the number of massive neutrino species and on the sum of their masses rriv from WMAP data 
alone, WMAP data and HST measurements of the Hubble constant, WMAP and SPT data, WMAP and SNLS3 data, and 
WMAP plus SPT plus SNLS3 data. 



here, the strong degeneracy existing among the two neu- 
trino parameters m„ and N c g makes mandatory the 
combination of several data sets. The constraints from 
the combination of WMAP and HST measurements be- 
come one order of magnitude weaker than in the previ- 
ous two scenarios (N c & fixed or rriy fixed), although if 
more data sets are added in the analysis, the situation 
improves significantly. Given that we know the sum of 
neutrino masses is non-zero, this means that constraining 
the number of neutrino species while neglecting the sum 
of their masses could lead to completely wrong results, 
and thus in general caution is advised. In the reverse sce- 
nario, i.e. constraining m„ while keeping N c s fixed, the 
same danger exists, but the major difference is that N e g 
may well be equal to its standard model value, while we 
know Yjra v is non-zero. 

To summarize the result with our maximal combina- 
tion of data sets, the joint constraints on the effective 
number of massive neutrino species and on the sum of 
their masses from CMB, supernovae and the new BAO 
data are 3. 71^0^2 an d J2 m v < 0.51 eV, respectively 
(both at the 95% CL). Near future cosmological data, 



as CMB measurements from the Planck mission, com- 
bined with other probes, are expected to settle the issue 
of whether there exists cosmological compelling evidence 
for extra radiation. 
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FIG. 3: Left panel: the yellow contours show the 68% and 95% CL constraints in the (N e g, ^ m v ) plane from our basic WMAP 
data set combined with old BAO data. The blue contours show the results from the combination of WMAP, SPT and HST 
measurements, while the red contours depict the results from the combination of WMAP, SPT and SNLS3 measurements. 
Finally, the green contour denotes the constraints from the combination of WMAP, SPT, SNLS3 and BAO data. Right panel: 
as in the left panel but using recent BAO2012 data. 



11 



[51] F. Beut ler et al, Mon. N ot. Roy. Astron. Soc. 416, 3017 
(2011) [arXiv:1106.3366l [astro-ph.CO]]. 



